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Thermal transport measurements have been performed on single-crystalhne Co-doped 
BaFe2As2 down to 0.1 K and under magnetic fields up to 7 T. Significant peak anoma- 
lies are observed in both thermal conductivity and thermal Hall conductivity below Tc as 
an indication of the enhancement of the quasiparticle mean-free path. Moreover, we find a 
sizable residual T-lincar term in thermal conductivity, possibly due to a finite quasiparticle 
density of states in the superconducting gap induced by impurity pair breaking. Our find- 
ings support a pairing symmetry compatible with the theoretically predicted sign-reversing 
s-wave state. 
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The symmetry of the order parameter is essential for identifying the superconducting 
pairing mechanism. In conventional superconductors (e.g., Al and Pb), the effective electron 
interaction is mediated by phonons, which gives rise to the isotropic s-wave pairing symmetry. 
On the other hand, electron pairs glued by magnetic interactions form unconventional pairing 
states: p-, d-wave states and so on. So far, such unconventional superconducting states have 
been found in a number of materials on the boarder between magnetism and superconductiv- 
ity.^ These findings suggest that a system close to magnetic instability provides a fertile field 
for unconventional superconductivity. A new family of superconductors containing layers of 
iron pnictides bear resemblance to unconventional superconductors such as high-Tc cuprates 
with a two-dimensional electronic structure and a magnetic order proximity to the supercon- 
ducting phase. ^'^ Therefore, an exotic superconducting pairing state can be naively expected 
in this system. In fact, an intriguing pairing state of sign-reversing s-wave symmetry has been 
theoretically proposed.^' ^ 

Here, we report the first thermal transport evidence of a novel pairing state in Co-doped 
BaFe2As2. In particular, we find a sizable residual T-linear term of the thermal conductivity, 
possibly due to the impurity-induced in-gap state. In addition, significant peak anomalies 
are observed in both thermal conductivity and thermal Hall conductivity originating from 
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the prominent enhancement of the quasiparticle (QP) mean-free path below T^- The field 
dependence of the delocalized QP density of states is apparently different from that of nodal 
gap excitation. These observations all point to the fully gapped sign-reversing s-wave state. 

Single-crystalline samples with the composition Ba(Feo.93Coo.o7)2As2 were grown by the 
FeAs/CoAs self-flux method.^ The substituted Co atoms donate extra electrons to FeAs lay- 
ers as itinerant carriers without creating localized moments.^ For thermal transport measure- 
ments, two different crystals with dimensions of 1 x 0.4 x 0.04 mm'^ were used. A one-heater- 
two-therniometcr steady-state method was used to measure thermal conductivity down to 0.1 
K and up to 7 T. Heat current was always aligned within the ah plane of the sample, and 
the magnetic field is applied along the c-axis. The samples were cooled in a magnetic field 
to ensure field homogeneity. We used Cernox and Ru02 thermometers above and below 1 K, 
respectively. The thermometers were thermalized on the sample by gold wires held by spot 
welding, providing a good thermal contact with a low electrical contact resistance, Rc < 50 
mil, at 300 K. In fact, we confirmed ohmic thermal response with the applied power down to 
the lowest temperature. The same contacts and gold wires were used to measure the resistivity 
of the sample by a standard four-contact method. 

First, we present the temperature dependence of the thermal conductivity Kxx{T)/T under 
several magnetic fields in Fig. 1. The arrow indicates Tc at zero field determined from the 
resistivity measurement shown in the upper inset of Fig. 4. Interestingly, Kxx(T)/T increases 
below Tc and reaches a maximum at 12 K (-^ Tc/2). In addition, the peak is suppressed by 
applying a magnetic field. In general, thermal conductivity is composed of the electronic term 

and the phononic term k?^x\ i^-xx = f^%x + '^Sx- L^t us estimate k%.j. by assuming that the 
Wicdemann-Franz law is valid at T^. We obtain ^.^./T = L^/p = 1.6 x 10"^ W/K^ni ~ 8.4 % 
of Kxx/T, indicating a predominant phononic contribution above Tc. Here, Lq = 2.44 x 10~^ 
Wr2/K^ is the Lorentz number and p = 150 p,Qcm. Interestingly, a similar increase and its 
field suppression are observed in hole-doped Bai_xKxFe2As2 as well as in unconventional 
superconductors such as YBa2Cu307_5^ and CeCoIns,^'' originating from the enhancement of 
the QP mean-free path in the superconducting state. It should be noted that the enhancement 
factor Kxx{Tc/2) / Kxx{Tc) in Bai_a;Ka,Fe2As2^ is two times larger than that in our result. The 
origin of this discrepancy will be discussed below. 

To examine the electronic contribution to heat transport, we quantitatively estimated 
and in the superconducting state from the field dependence of the thermal conductivity 
i^xxiB) (inset of Fig. 1). The data can be well described by the vortex-scattering model^'^^ 

'^xx(r,5) = ^^^|^ + /.P^(T), (1) 

where a(T) is an inverse field scale. In the analysis, /ix^(T) and k%^{T) are fixed to be 
field— independent by assuming that the condensate amplitude is nearly unaffected by mag- 
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Fig. 1. (Color online) Temperature dependence of thermal conductivity Kxx{T)/T under several 
magnetic fields for BJ^q and q\\ab. Inset: Field dependence of normalized thermal conductivity 
Kxx{B) / Kxx{0) at fixed temperatures. Solid lines are fittings by the vortex scattering model (see 
text). 

netic field at B < 7 T}^ Figure 2(a) shows the temperature dependences of k%^{T)/T and 
K^{T)/T determined by fitting the data using eq. (1). It is noteworthy that k%^{T) /T exhibits 
a maximum at around Tc/2 similar to Kxx{T)/T. This implies that the peak in Kxx/T at zero 
field is associated with k%x- The fittings also provide a(T), which is proportional to the QP 
mean-free path I through the relation a{T) = Zo"tr/^o/^ where the cross section cxtr is roughly 
equal to the coherence length ^ and $0 is the flux quantum. In the inset of Fig. 2(b) (left 
axis), we show the temperature dependence of a(T), which exhibits a steep increase below 
Tc- This reflects an enhancement of I in the superconducting state. At 2 K, / is found to be ~ 
1100 A with ^ ~ 34 A from H^^^ ~ 30 T.^^ 

Alternatively, the enhancement of the QP mean-free path can be confirmed from the ther- 
mal Hall conductivity K,xyi which is a powerful probe for QPs since it is purely electronic. Fig- 
ure 3(a) presents the field dependence of |K^y|(i?) at fixed temperatures. The sign of Kxy{B) is 
negative, consistent with that in electron doping opposite to the hole-doped Bai-xKxFe2As2.^ 
In Fig. 2(b), we plot the temperature dependence of the initial slope \Kxy\/B= lim \Kxy\/B. 
As is clearly seen, \Kxy\{T) / B peaks at around Tc/2 similar to the temperature dependence 
of KxxiT) /T ■ This directly indicates the enhancement of the electronic contribution to heat 
transport below Tc. In addition, the QP mean-free path is provided by the thermal Hall angle 
tan@ = l^xyl/K'^x because tanQ/B oc / in the weak field limit. ^'^ We plot the temperature 
dependence of tanQ/B in the inset of Fig. 2(b) (right axis). Clearly, tanO/S increases with 
decreasing temperature, as we found in a{T) cx /. Therefore, all data consistently point to 
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Fig. 2. (Color online) (a) Temperature dependence of electronic term k%^{T)/T and phononic term 
K^{T)/T of thermal conductivity, (b) Temperature dependence of initial slope of thermal Hall 
conductivity \Hp^y\/B= lim \Hixy\/B. Inset: Temperature dependence of a(T) (left axis) and ther- 
mal Hall angle tanO/B (right axis). 



the striking enhancement of / below T^. This, in turn, suggests that I in the normal state is 
suppressed by the inelastic scattering possibly due to antiferromagnetic (AF) fluctuations, as 
discussed in the microwave conductivity measurements.^^ On the other hand, in the super- 
conducting state, electrons condense into Cooper pairs and their number decreases rapidly 
below Tc. This in turn gives rise to a reduction in the scattering cross section of QPs, and 
hence / increases below Tc. The presence of AF fluctuations has also been indicated by recent 
NMR measurements.^ 

The thermal Hall conductivity Kxy also provides the density of states of delocalized QPs, 
^dc\{E), using the conjectures oc N^^iiE)! and \Kxy\/ k'^^B oc /. A plot of k^,^.'^ B /\Kxy\ as 
a function of field reveals the field dependence of A'dei (E) . Note that the precise estimation of 
Niic\{E) from specific heat is rather difficult in iron pnictides because of the contribution of the 
nuclear Schottky anomaly. As seen in Fig. 3(b), K%^'^B/\Kxy\ shows a weak field dependence 
within the experimental error of \Kxy\, in contrast to the strong field dependence observed 
in nodal superconductors.^'' One may expect such a weak field dependence in fully gapped 
superconductors because QPs localized around vortex cores do not contribute to thermal 
conductivity at low fields B ^ Hc2- 

Next, we discuss the low-temperature thermal conductivity to gain insight into the su- 
perconducting pairing symmetry. We present the temperature dependence of the thermal 
conductivity Kxx{T) down to 0.1 K under zero magnetic field in Fig. 4. The lower inset of 
Fig. 4 shows the low-temperature part of the Kxx/T vs plot. The straight line is a fit to 
Hxx/T = kq/T + bT'^, where kq/T is the residual T-linear term extrapolated to T — > K. The 
best fit was obtained with kq/T = 1.2 x 10~^ W/K^m and b = 0.33 W/K^m. What is surpris- 
ing is that the residual T-linear term kq/T amounts to as much as half of the normal— state 
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Fig. 3. (Color online) Field dependence of (a) thermal Hall conductivity \Kxy\{B) for BJ^q and q\\ab 
and (b) k^^'^ B /\Kxy\ which is proportional to the density of states of delocalized quasiparticles at 
fixed temperatures. 



thermal conductivity Kn/T, which is estimated from the Wiedemann-Franz law as Hn/T = 
Lq/po = 1.9 X 10"'^ W/K-^m. Here, po = 130 /iJlcm is the residual resistivity at T = K 
obtained by assuming that p decreases linearly against temperature below Tc (see dashed line 
in the upper inset of Fig. 4). On the other hand, using the mean acoustic phonon velocity 
(vg) = 2400 m/s and the phonon specific heat coefficient /? = 8.98 J/K^m^ obtained from 
the parent compound of BaFe2As2,^^ the slope of 6 = |/3(us)/ph yields /ph = 46 pm, which 
is the same order of magnitude as the smallest crystal dimension, namely, 40 pm thickness 
of the sample. This implies that the low-temperature thermal conductivity is dominated by 
phonons. 

In the superconducting state, the presence of the residual T-linear term kq/T can be at- 
tributed to the impurity bound states formed by the interference of particle- and hole-like 
excitations, which undergo Andreev scattering, evoking sign changes of the order parameter 
as a result of unconventional pairing and impurity scattering. Moreover, in the nodal super- 
conductors, kq/T takes a universal value independent of impurity concentration because the 
quasiparticles are both generated and scattered by the impurities. In fact, the impurity con- 
centration independent kq/T has been observed in high-Tc superconductors,^^ although its 
universality is still under debate. Theoretically, the universal thermal conductivity kqq/T 
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Fig. 4. (Color online) Temperature dependence of thermal conductivity Kxx{T) for q\\ab under zero 
field. Upper inset: Temperature dependence of resistivity p{T) under zero field. Lower inset: k^x/T 
vs plot under zero field. The solid line represents a fit to the data by Kxx/T — kq/T + bT^. 



for the nodal superconductor is explicitly expressed as 

-T^ = — NfVf X — , (2) 

where Nj is the normal density of states, Vf is the Fermi velocity, Aq is the maximum am- 
plitude of the gap, fi is the slope of the gap at the node on the circular Fermi surface, and 
a is ^ for the 2D order parameter with lines of nodes. Using the Wiedemann-Franz law, 
eq. (2) can be written as ^ = (f ^^)^-^'' By assuming a d-wave gap structure, we obtain 
kqo/T = 0.4 X 10-3 W/K^m for Ba(Feo.93Coo.o7)2As2 using fx = 2 with ^ = 34 A and Z = 
1100 A. The determined value is one order of magnitude smaller than the experimental value 
of hq/T. Although we cannot exclude the possibility of the nodal superconducting state on 
the basis of only kqq/T, given the leading phononic contribution in low-temperature k^x, as 
evidenced by the term and the weak field dependence of N^^^d, nodal superconductivity is 
highly unlikely in Ba(Feo.93Coo.o7)2As2. Here, a question is raised as to how one can take a 
finite T-linear term in f^xx 

/T without the nodal pairing state. 
In addition to the nodal superconducting state, one can also expect a residual T-linear 
term in a novel s-wave pairing state. At first glance, the residual T-linear term seems to be 
incompatible with the s-wave state on the basis of Anderson's theorem. On the contrary, it 
is possible to realize the finite value via impurity interband scattering between different Fermi 
surfaces with opposite signs of order parameters with full gaps. In fact, it is theoretically 
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proposed that the AF fluctuations with Q ~ (^^ 0) resulting from the interband nesting 
between the hole and electron pockets realize a fully gapped sign-reversing s-wave state {s±- 
wave state) in the iron pnictide superconductors.^'^ In this state, the impurity- induced density 
of states (DOS) in the gaps is predicted to emerge when interband scattering is enhanced by 
the introduction of impurities. However, it should be emphasized that such an in-gap state 
can be expected only if a strong impurity pair breaking involving a large reduction in T^, 
which amounts to 10 K, occurs. 

Here, we discuss whether such a large reduction in Tc occurs in Co-doped BaFe2As2 to 
examine the possibility of the s-j--wave state in this system. Given that Tc = 25 K of Co-doped 
BaFe2As2 is about 10 K lower than that of the K-doped one, Tc = 37 K,^ the difference in Tc 
can be explained by the difference in the strength of the pair breaking caused by the dopant. 
This is because the Co atoms substituted into the conducting layer can act as stronger pair 
breakers than the K atoms doped into the block layer. In fact, the enhancement of the thermal 
conductivity HxxiTcf^) / f^xxiTc) for Ba(Feo.93Coo.o7)2As2 is reduced to half of that for K-doped 
BaFe2As2.^ In addition, the resistivity at Tc for Ba(Feo.93Coo.o7)2As2 {p = 150 ^uilcm) is three 
times larger than that for the K-doped one (p = 50 /iilcm).^ These results indicate that the 
scattering rate of Ba(Feo.93Coo.o7)2As2 is practically enhanced by Co doping. In this case, one 
may expect the absence of an in-gap state or a smaller in-gap state in K-doped BaFe2As2 with 
weak scattering by the K atoms. Recently, a negligible residual T- linear term in Kxx/T was 
found in Bai_2,K2,Fe2As2.^^ Moreover, the spin-lattice relaxation rate 1/Ti of the Co-doped 
BaFe2As2 appears to level off at low temperatures,^ while that of the K-doped one varies 
close to T^ all the way down to the lowest measured temperature.^'^ These results further 
explain the role of Co atoms as strong pair breakers and support our scenario. Thus, our 
findings all point to the fully gapped s±-wavc state in Ba(Feo.93Coo.o7)2As2, being consistent 
with the theoretical prediction as well as the resonance peak observed by the inelastic neutron 
scattering measurement.^^ 

In summary, we uncovered an intriguing pairing state of the fully gapped sign-reversing 
s-wave state in the Co-doped BaFe2As2 by thermal transport measurements. Further ex- 
periments on samples with different Co doping levels are in progress to clarify the effect of 
inhomogeneity in the sample and the crucial role of non-magnetic impurities to the residual 
in-gap state. 
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